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Abstract. In this paper we report new observations and derive chemical abundances for a sample of 26 planetary 
nebulae (PN) located in the anticenter direction. Most of these nebulae are far away objects, located at galac- 
tocentric distances greater than about 8 kpc, so that they are particularly useful for the determination of the 
radial gradients at large distances from the galactic center. A comparison of the present results with previously 
determined abundances suggests that the radial abundance gradients flatten out at distances larger than about 
10 kpc from the center. 
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1. Introduction 

The presence of radial abundance gradients for several el- 
ement ratios in the Milky Way is now well established (see 
Maciel 1997, 2000 for a revision). Gradients can be derived 
from different kinds of objects, such as 0,B stars, HII re- 
gions or planetary nebulae, and seem to be a common 
property of spiral galaxies. In fact, there are reports of gra- 
dients in many other galaxies apart from the Milky Way 
from the analysis of HII regions fSearle H9711 Smith 1*11)751 
Blair et al. ITlfel Kennicutt & Garnett IT55rjl Ferguson et 
aLHSSJ. 

The study of radial gradients involves one of the most 
important constraints for galactic evolution models, as 
they can provide at least three crucial informations: i) 
the magnitude of the gradient itself, generally measured 
in dex/kpc; ii) possible variations in the gradient for dif- 
ferent galactocentric distances; iii) the time variation in 
the gradients. 

Regarding the space variations of the gradients, 
presently available results are somewhat controversial. 
Some results for the Milky Way based on stars (Kaufer 
et al. ITSM*) and HII regions (Esteban & Peimbert ITM>1 
Vilchez & Esteban Il99rjf) suggest decreasing gradients for 
galactocentric distances R > 6 kpc. This is confirmed by 
a recent analysis based on planetary nebulae (Maciel & 
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Quireza ,1999), and is also supported by chemical evo- 
lution models by Chiappini et al. i 1 997; and Matteucci 
(2000). However, no significant variations in the gradient 
are detected from a sample of O, B stars by Smartt |2000 
and references therein) , and a flattened gradient near the 
solar neighbourhood has been suggested by Andrievsky et 
al. 1)20021 and references therein) on the basis of a sample 
of cepheid variables. 

The shapes of the gradients in spiral galaxies are also 
subject to discussion. Some flattening in the outer regions 
has been proposed (Diaz [1989), while some recent results 
are consistent with an exponential gradient with a con- 
stant slope (Henry & Howard 11995). 

Planetary nebulae can play an important role in this 
context (Maciel & KoppenUHMl Peimbert & Carigi UMSjl . 
as nebular abundances can be accurately derived, even 
for objects at large heliocentric distances, which is par- 
ticularly useful for objects in the direction of the galactic 
anticenter (Maciel & Quireza 1999). Also, these objects 
can be classified according to the masses and ages of their 
progenitors, so that they are also useful in the study of 
the time variations of the gradients (Maciel & Costa[2003 
Maciel et al. l2003fl . The derivation of chemical abundances 
for a homogeneous and wide sample of PN at large galac- 
tocentric distances, located basically outside the solar cir- 
cle of about 8 kpc, would therefore provide the tools for 
the assessment of the gradient properties at these large 
distances. 
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In this work, we report the chemical abundances for 
a sample of PN towards the galactic anticenter, where 
most of the objects placed at large galactocentric distances 
are located. They were selected from the ESO-Strasbourg 
Galactic PN Catalog (Acker et al. I1992[) . and the sample 
consists of 26 objects, all of them observed, reduced and 
analysed homogeneously. 



2. Observations 

The observations were carried out in two different ob- 
servatories: at ESO/La Silla, and at LNA/Brazil. At 
ESO we used the 1.52m telescope and the Boiler & 
Chivens Cassegrain spectrograph, in two observational 
runs: dec/1999 and dec/2000. We used grating #23 and 
CCD #38, which provide a dispersion of 2.2 A/pix, and 
spectral coverage of about 5000 A. The effective wave- 
length coverage was reduced to about 3500 A starting at 
4200 A, due to a low S/N ratio shortward of this wave- 
length. At LNA, in dec/1995, we used the 1.60m telescope, 
a Boiler & Chivens Cassegrain spectrograph and a com- 
bination grating+CCD, which gives a dispersion of 4.5 
A/pixel and spectral coverage of about 4500 A. An east- 
west long slit of 2 arcsecs width was used for all oberva- 
tions. Table 1 gives the logbook of the observations and 
also the adopted distances d (kpc) to the Sun and galacto- 
centric distances R (kpc), which we will use in Section 4. 
Flux calibration was secured through the observation of 
at least three spectrophotometric standards in each night. 
Each object was observed at least twice, and line fluxes 
were derived from the average of all good measures for 
each object. 

Data reduction was performed using the IRAF pack- 
age, and followed the standard procedure for long slit spec- 
troscopy: correction of bias, dark and flat-field, extraction 
of the spectral profile, wavelength calibration and flux cal- 
ibration. Atmospheric extinction was corrected through 
the mean coefficients for the observatory and season. Line 
fluxes were calculated assuming gaussian profiles, and a 
gaussian de-blending routine was used when necessary. 
Fluxes calculated for all diagnostic lines are listed in 
Table 2 (available electronically at the CDS 1 ). Interstellar 
reddening was derived from the difference between mea- 
sured and predicted values for the Balmer ratio Ha/H/3, 
assuming Case B (Osterbrock 1989 ), and adopting the in- 
terstellar reddening law by Cardelli et al. H1989[) . with 
Ry=3.1. The derived E(B-V) values are given in Table 3 
for each nebula. 

A preliminary analysis of the objects observed in 1995 
was already published (Costa et al. 1997); however, for the 
present work we rederived the physical parameters and 
abundances for the whole sample, in order to assure the 
homogeneity of the results. 



3. Physical parameters and abundances 

Electron densities n e were derived from the 
[SII]A6717/6731 ratio, while electron temperatures 
T e were estimated from both the [OIII]A4364/5007 and 
[Nil] A 5755/6584 line ratios. When both line ratios were 
available, the [OIII] temperature was used to derive 
the ionic concentrarions of higher ionization potential 
ions such as 0+ 2 , Ar+ 2 , Ar+ 3 , Ne+ 2 , and the [Nil] 
temperature for lower ionization potential ions like + , 
N + , S + . In the cases when only one temperature was 
available, it was used for all ions. For the objects A20 and 
NGC 1514 no temperatures were available, so that an 
average value of 10 4 K was adopted. Concerning densities, 
for these nebulae the [SII] ratio was not available, and a 
typical value of 5000 cm -3 was adopted. Table 3 shows 
the adopted electron densities (cm -3 ) and temperatures 
(K) for all objects. 

Ionic abundances were calculated using a three- 
level atom model by solving the statistical equilibrium 
equations, including radiative and collisional transitions. 
Elemental abundances were then derived using ionization 
correction factors to account for unobserved ions of each 
element. The atomic data and ICFs used in this work are 
same as adopted by Escudero & Costa (2 001fl . to which 
the reader is referred for details. The adopted formulae for 
the ICFs are given below. 
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For the cases where [SIII] 6312 is not available, we adopted 
the expression given by Kingsburgh & Barlow (1993) to 
estimate the S ++ abundance. As in Escudero & Costa 
(2001), we calculated helium abundances using the re- 
combination theory, with recombination coefficients from 
Pcquignot et al. ( 1991 ), and correcting the Hel abundance 
from collisional effects using the correction terms from 
Kingdon & Ferland l|1995fl . He I abundances were deter- 
mined as weighted averages from the 4471, 5876 and 6678 
lines, using the line fluxes as weights. He II abundances 
were derived from the 4686 line. The obtained ionic abun- 
dances are given in Table 4, where "E-05" means 10~ 05 , 
etc. 
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Table 1. Log of the observations and distances 



Object 


PN G 


RA 


DEC 


Place 


& Date 


(/ 


R 


A 12 


198.6-06.3 


06 


02 


23 


09 39 03 


ESO - 


Dec. 


1999 


2.0 


9.5 


A 18 


216.0-00.2 


06 


56 


15 


-02 53 12 


ESO - 


Dec. 


1999 


1.6 


8.9 


A20 


214.9+07.8 


07 


23 


01 


01 45 37 


ESO - 


Dec. 


2000 


2.0 


9.3 


H 3-75 


193.6-09.5 


05 


40 


44 


12 21 16 


ESO - 


Dec. 


1999 


2.4 


9.9 


IC 2003 


161.2-14.8 


03 


56 


22 


33 52 30 


LNA - 


Dec. 


1995 


4.7 


12.0 


IC 2165 


221.3-12.3 


06 


21 


43 


-12 59 10 


LNA - 


Dec. 


1995 


2.0 


9.2 


J 320 


190.3-17.7 


05 


05 


34 


10 42 21 


LNA - 


Dec. 


1995 


6.1 


13.4 


J 900 


194.2+02.5 


06 


25 


57 


17 47 27 


LNA - 


Dec. 


1995 


2.8 


10.3 


K 1-7 


197.2-14.2 


05 


31 


48 


06 56 09 


LNA - 


Dec. 


1995 


3.8 


11.2 


K 2-1 


173.7-05.8 


05 


03 


09 


30 50 03 


ESO - 


Dec. 


2000 


1.1 


8.7 


K 3-66 


167.4-09.1 


04 


36 


37 


33 39 30 


LNA - 


Dec. 


1995 


6.7 


14.1 


K 3-70 


184.6+00.6 


05 


58 


45 


25 18 44 


ESO - 


Dec. 


1999 






M 1-6 


211.2-03.5 


06 


35 


45 


-00 05 41 


ESO - 


Dec. 


1999 


2.6 


9.9 


M 1-7 


189.8+07.7 


06 


37 


21 


24 00 30 


ESO - 


Dec. 


1999 


5.9 


13.4 


M 1-8 


210.3+01.9 


06 


53 


34 


03 08 24 


LNA - 


Dec. 


1995 


3.4 


10.7 


M 1-9 


212.0+04.3 


07 


05 


19 


02 46 59 


LNA - 


Dec. 


1995 


4.9 


12.0 


M 1-13 


232.4-01.8 


07 


21 


15 


-18 08 36 


LNA - 


Dec. 


1995 


5.3 


11.6 


M 1-14 


234.9-01.4 


07 


27 


56 


-20 13 23 


LNA - 


Dec. 


1995 


4.0 


10.4 


M 1-16 


226.7+05.6 


07 


37 


19 


-09 38 48 


LNA - 


Dec. 


1995 


5.5 


12.0 


M 1-17 


228.8+05.3 


07 


40 


22 


-11 32 30 


LNA - 


Dec. 


1995 


7.4 


13.6 


M 1-18 


231.4+04.3 


07 


42 


06 


-14 21 41 


ESO - 


Dec. 


1999 


4.4 


10.9 


M 3-2 


240.3-07.6 


07 


14 


50 


-27 50 20 


ESO - 


Dec. 


1999 


8.8 


14.1 


M 3-3 


221.7+05.3 


07 


26 


35 


-05 21 57 


ESO - 


Dec. 


1999 


5.7 


12.4 


NGC 1514 


165.5-15.2 


04 


09 


22 


30 46 00 


ESO - 


Dec. 


2000 


0.8 


8.3 


SA 2-21 


238.9+07.3 


08 


08 


43 


-19 13 59 


ESO - 


Dec. 


1999 


3.3 


9.7 


YC 2-5 


240.3+07.0 


08 


10 


42 


-20 31 33 


ESO - 


Dec. 


1999 


7.7 


13.2 



Table 2. Line fluxes (available electronically) 



Ion 


A 12 


A 18 


A 20 


H3-75 


IC 2003 


IC 2165 


J 320 


J 900 


H 7 4340 


54.0 


34.0 


27.0 


68.0 


26.0 


53.0 


63.0 


59.0 


[O III] 4363 


8.7 


10.6 


3.8 


13.4 


11.5 


14.3 


13.6 


9.6 



3.1. Abundances and distance-independent correlations 

Derived elemental abundances are given in table 5. 
The helium abundance (He/H) is relative to hydrogen, 
and the remaining elemental abundances are given as 
e(X)=log(X/H)+12. The N/O ratio is also given, as this 
ratio is useful in order to classify the nebulae, and is also 
related to the progenitor star mass. We have also exam- 
ined the classification of the objects in our sample ac- 
cording to the classification scheme proposed by Peimbert 
(1978), which relies basically on elemental abundances. 
All nebulae in the sample are disk objects, belonging to 
types I, II and III. Subtypes Ha and lib (Faundez-Abans 
& Maciel lT987[) have also been included. Bulge (type V) or 
halo (type IV) objects were avoided (see the last column 
of Table 5). 

Average helium abundances in our sample are He/H = 
0.160 for type I PN and He/H = 0.107 for type II objects. 
For oxygen, we obtain e(0) = 8.42 and e(O) = 8.66, re- 



spectively, thus showing some indication of ON cycling in 
type I PN, presumably those with more massive progeni- 
tors in which the second dredge up process has operated 
(see for example de Freitas Pacheco 119931 Kingsburgh & 
Barlow IT^jl . 

The relative abundances can be better examined us- 
ing distance- independent correlations, as can be seen in 
Figure n which shows log(N/0)x He/H. Our new data 
are presented as solid dots, and the crosses show the re- 
sults of our previous work as given by Costa et al. (^996 ) , 
de Freitas Pacheco et al. (HHBHl HHHII CHH2I) , and Maciel et 
al. (1990). As can be seen in the previous papers, all ob- 
jects have been observed and analyzed in a homogeneous 
way, so that the total sample includes about 80 nebulae 
in the disk alone, forming one of the largest database of 
accurate abundances of planetary nebulae available. This 
is important, since distance independent correlations as 
shown by Fig. ^ are much better defined when a large 
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Table 3. Physical parameters 



Object 


E(B - V) 


T e [OIII] 


T e [Nil] 


7le[SII] 


A 12 


0.440 


11912 


9144 


255 


A 18 


0.969 


15736 


14317 


364 


A 20 


0.0 






5000 


H 3-75 


0.290 


12813 


9013 


119 


IC 2003 


0.301 


12080 


13276 


5144 


IC 2165 


0.417 


12519 


14075 


5715 


J 320 


0.248 


11944 


12815 


4915 


J 900 


0.558 


11320 


12513 


4016 


K 1-7 


0.144 


16123 


11329 


92 


K 2-1 


0.433 


11375 




493 


K 3-66 


0.660 


9057 


13143 


4725 


K 3-70 


1.530 


17500 


11643 


2361 


M 1-6 


1.260 


10558 


11306 


8426 


M 1-7 


0.129 


10244 


11016 


1135 


M 1-8 


0.374 


15563 


10307 


482 


M 1-9 


0.540 


10922 


10435 


4516 


M 1-13 


0.208 


8196 


9601 


994 


M 1-14 


0.757 


12805 


9987 


6784 


M 1-16 


0.450 


10657 


10351 


2311 


M 1-17 


0.650 


10008 


10525 


6030 


M 1-18 


0.088 


11224 


10896 


69 


M 3-2 


0.317 




10164 


230 


M 3-3 


0.250 


13826 


9867 


349 


NGC 1514 


0.410 






5000 


SA 2-21 


0.200 


12720 


15152 


351 


YC 2-5 


0.290 


11477 


11477 


5000 



sample is considered. It can be seen that our new data fit 
rather nicely in the observed correlation. 

The plot shows the expected trend of an increase in 
N/O as He/H increases, reflecting the enrichment of He 
and N due to dredge-up episodes, as predicted by stellar 
evolution models. Some N/O enrichment can also be ob- 
served by comparing the data in Fig. ^ ( or Table 5) with 
the corresponding values of the N/O ratio for HII regions 
of similar metallicity, measured by the O/H ratio, as given 
for example by Henry & Worthey (1999). 

Figures[3]and[3]show the correlations with the O/H ra- 
tio of the Ar/H and S/H abundances, respectively. These 
correlations are also well defined, and reflect the mass dis- 
tribution of PN progenitors, as younger, more massive ob- 
jects are expected to form from an interstellar medium 
richer in alpha-elements. The Ar abundance of the neb- 
ula YC 2-5 (Sa 2-22) shows a larger deviation from the 
trend as compared to the remaining objects, as indicated 
in Figure |2 This object shows a normal He abundance, 
and its oxygen abundance can also be considered as nor- 
mal, albeit larger than the average. If that is true, the 
argon abundance is probably a lower limit, which may be 
due to the fact that its electron density is not accurately 
determined. 



3.2. Comparison with previous results from the 
literature 

Most of the nebulae studied in this work are objects lo- 
cated at relatively large distances from the sun, for which 
accurate elemental abundances are being determined for 
the first time. Some objects, however, have been studied 
by other groups, so that it is interesting to compare our 
new results with previously determined abundances in the 
literature. 

Considering first the earlier work, the objects IC 2003, 
IC 2165, J320 and J900 are part of the detailed study 
by Aller & Czyzak l|1983[) based on image-tube scanner 
measurements and making use of theoretical models to 
represent individual line intensities. Regarding the Hc/H 
abundances, our data agree with their results within 4% 
for IC 2003, J320 and J900, with a larger discrepancy of 
about 38% for IC 2165, for which the abundance by Aller 
& Czyzak ( 1983 ) is higher than ours. Our electron temper- 
atures for this object do not differ from the adopted values 
by Aller & Czyzak (I983JI and more recently by Kholtygin 
( 199SJ) by more than 8%. However, a larger discrepancy in 
the electron temperatures of 11 to 16% and in the helium 
abundance again of 38% is apparent for this nebula from 
the detailed work by Hyung ( 1994), so that our ionization 
correction factor may not account completely for the to- 
tal abundance, and our result is probably underestimated. 
Considering now the heavy elements O, N, S and Ar, for 
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Table 4. Ionic Abundances 



Ubject 


He I 


He II 




N II 




Q TT 
b 11 




S III 




O II 




C\ TTT 

U 111 




Ar III 


A 12 





111 





008 


8 


45E- 


05 


2 


82E-06 


3 


23E- 


06 


6 


21E- 


04 


1 


69E-04 


9 


03E- 


07 


A 18 





141 





010 


4 


27E- 


05 


5 


25E-07 


1 


87E- 


06 


4 


43E- 


05 


4 


60E-05 


5 


25E- 


07 


A 20 





042 





083 








1 


59E-07 


7 


68E- 


06 


4 


63E- 


05 


1 


66E-04 


9 


93E- 


08 


H 3-75 





069 





003 


2 


50E- 


05 


9 


15E-07 


2 


15E- 


06 


1 


45E- 


04 


1 


77E-04 


6 


68E- 


07 


IC 2003 





039 





054 


1 


50E- 


06 


7 


99E-08 


1 


25E- 


06 


1 


07E- 


05 


2 


04E-04 


4 


63E- 


07 


IC 2165 





047 





029 


2 


51E- 


06 


9 


38E-08 


1 


75E- 


06 


9 


28E- 


06 


2 


09E-04 


4 


98E- 


07 


J 320 





101 





006 


2 


25E- 


07 


2 


33E-08 


9 


80E- 


07 


4 


23E- 


06 


2 


57E-04 


4 


06E- 


07 


J 900 





056 





038 


4 


22E- 


06 


1 


44E-07 


1 


33E- 


06 


2 


96E- 


05 


2 


47E-04 


4 


75E- 


07 


K 1-7 





100 





016 


2 


79E- 


05 


1 


07E-06 


1 


11E- 


06 


1 


16E- 


04 


8 


60E-05 


5 


33E- 


07 


K 2-1 





056 





066 


1 


89E- 


05 


3 


25E-07 


2 


41E- 


06 


1 


63E- 


04 


2 


52E-04 


1 


02E- 


06 


K 3-66 





083 






9 


26E- 


06 


8 


93E-08 


2 


03E- 


06 


4 


83E- 


05 


1 


21E-04 


6 


86E- 


07 


K 3-70 





138 





023 


1 


16E- 


04 


9 


61E-07 


2 


14E- 


06 


1 


05E- 


04 


7 


30E-05 


7 


09E- 


07 


M 1-6 





081 






2 


42E- 


05 


1 


91E-07 


1 


73E- 


06 


1 


88E- 


04 


3 


12E-05 


4 


34E- 


07 


M 1-7 





132 





005 


5 


97E- 


05 


9 


79E-07 


5 


58E- 


06 


2 


26E- 


04 


3 


85E-04 


3 


45E- 


06 


M 1-8 





122 





034 


1 


02E- 


04 


1 


20E-06 


1 


27E- 


06 


2 


56E- 


04 


1 


28E-04 


1 


10E- 


06 


M 1-9 





062 





001 


1 


12E- 


05 


1 


36E-07 


1 


07E- 


06 


8 


93E- 


05 


1 


18E-04 


4 


17E- 


07 


M 1-13 





122 





012 


1 


32E- 


04 


1 


23E-06 


5 


56E- 


06 


4 


80E- 


04 


7 


79E-04 


4 


74E- 


06 


M 1-14 





086 






1 


59E- 


05 


2 


23E-07 


5 


95E- 


07 


1 


53E- 


04 


4 


79E-05 


3 


64E- 


07 


M 1-16 


o 


104 


o 


023 


1 


40E- 


04 


3 


42E-07 


1 


18E- 


06 


2 


40E- 


04 


3 


46E-04 


1 


45E- 


06 


M 1-17 





069 





Oil 


3 


86E- 


05 


1 


50E-06 


6 


29E- 


06 


1 


61E- 


04 


4 


67E-04 


1 


57E- 


06 


M 1-18 





136 





017 


9 


21E- 


05 


1 


14E-06 


5 


97E- 


06 


1 


75E- 


04 


1 


91E-04 


1 


38E- 


06 


M 3-2 





134 





077 


1 


80E- 


04 


1 


21E-06 


6 


03E- 


07 


8 


95E- 


05 


1 


74E-05 


2 


51E- 


07 


M 3-3 





100 





029 


1 


79E- 


04 


3 


00E-07 


1 


22E- 


04 


1 


51E- 


04 


1 


22E-04 


8 


38E- 


07 


NGC 1514 





076 





005 








1 


63E-07 


3 


40E- 


06 


1 


96E- 


04 


4 


15E-04 


4 


48E- 


07 


SA 2-21 





110 





019 


1 


67E- 


05 


5 


72E-07 


2 


38E- 


06 


2 


99E- 


05 


1 


44E-04 


9 


65E- 


07 


YC 2-5 





017 





054 
























2 


66E-04 


1 


74E- 


07 



o 



Lf) 

d 



o 



LD 



o 



D.05 



+ 



_±_ 



0.' 



0.15 



He/H 



M3-2 



0.2 



Fig. 1. log(N/0) as a function of the helium abundance. Our new results are presented as solid dots, while the crosses 
show the abundances derived in our previous work, as quoted in the text. 



IC 2003, IC 2165 and J900 the agreement is good, with 
an average discrepancy of a factor 1.3. The worst case 
is the Ar abundance of IC 2165, reaching about a factor 
1.8, which is probably a result of the weaker lines of this 
element, or may be due the the adopted electron temper- 



ature. For these elements the difference of our results to 
the analysis by Hyung l|1994f) is of a factor 1.4 in average. 
Note however that our oxygen abundance agrees very well 
with the more recent work of Kholtygin ( 19SSJ, as we will 
see later on this section. For J320 the agreement is very 
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Table 5. Abundances 



Object . He/H e(O) e(N) e(S) e(Ar) log(N/0) Type 



A 12 


0.119 


8.93 


8.06 


6.78 


6.78 


-0.87 


Ila 


A 18 


0.152 


7.99 


7.97 


6.40 


6.17 


-0.02 


I 


A 20 


0.125 


8.80 




6.66 


6.86 




Ila 


H 3-75 


0.071 


8.52 


7.76 


6.52 


6.23 


-0.76 


lib 


IC 2003 


0.094 


8.71 


7.86 


6.53 


6.19 


-0.85 


lib 


IC 2165 


0.076 


8.55 


7.98 


6.64 


6.05 


-0.57 


lib 


J 320 


0.107 


8.44 


7.17 


6.45 


5.77 


-1.27 


lib 


J 900 


0.094 


8.67 


7.82 


6.42 


6.08 


-0.85 


lib 


K 1-7 


0.115 


8.37 


7.75 


6.36 


6.29 


-0.62 


lib 


K 2-1 


0.122 


9.19 


8.26 


6.62 


6.28 


-0.93 


Ila 


K 3-66 


0.083 


8.23 


7.51 


6.39 


6.11 


-0.72 


lib 


K 3-70 


0.161 


8.32 


8.36 


6.51 


6.43 


0.04 


I 


M 1-6 


0.081 


8.34 


7.45 


6.28 


6.61 


-0.89 


lib 


M 1-7 


0.137 


8.80 


8.23 


6.86 


6.88 


-0.57 


Ila 


M 1-8 


0.156 


8.69 


8.29 


6.41 


6.75 


-0.40 


Ila 


M 1-9 


0.063 


8.32 


7.42 


6.11 


6.00 


-0.90 


lib 


M 1-13 


0.134 


9.14 


8.58 


6.88 


7.05 


-0.56 


Ila 


M 1-14 


0.086 


8.30 


7.32 


5.91 


6.31 


-0.98 


III 


M 1-16 


0.127 


8.86 


8.62 


6.23 


6.60 


-0.24 


Ila 


M 1-17 


0.081 


8.86 


8.24 


6.98 


6.52 


-0.62 


lib 


M 1-18 


0.153 


8.62 


8.34 


6.88 


6.60 


-0.28 


I 


M 3-2 


0.212 


8.23 


8.53 


6.27 


6.51 


0.30 


I 


M 3-3 


0.129 


8.54 


8.62 


6.49 


6.51 


0.08 


I 


NGC 1514 


0.081 


8.81 




6.61 


5.97 




lib 


SA 2-21 


0.129 


8.31 


8.05 


6.61 


6.26 


-0.22 


Ila 
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9.05 
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Fig. 2. Argon abundances as a function of the O/H ratio. The symbols are as in Fig. 1. 



good for oxygen, which is well measured in both works, 
and poorer for the other elements. It should be noted that 
the sulfur and nitrogen abundances in the work of Aller & 
Czyzak (1983 ) are of quality C or D, and that some recent 
work that we will quote later also produces discrepant re- 



sults, which may be an indication of some variation inside 
the nebula. 

In the spectrophotometry surveys performed by 
Koppen et al. !|1991f> and Cuisinier et al. IjlHHfijl . about 140 
planetary nebulae have been studied on the basis of good 
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.5 9 
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Fig. 3. Sulfur abundances as a function of the O/H ratio. The symbols are as in Fig. 1. 



quality spectra in which plasma diagnostics have been 
tested with theoretical photoionization models. The first 
survey has four objects in common with our work, namely 
J320, J900, Ml-6 and M3-2. The He/H abundances for 
these objects are in agreement with our results within an 
average of 26%. The worst cases are J320 and Ml-6. In 
the former, the helium abundance was probably under- 
estimated by Koppen et al. l|1991f) . as our value agrees 
very well with Aller & Czyzak Q1983[l and with the re- 
cent work of Milingo et al. 1|2002[1 , as we will see later on 
in this section. For Ml-6, the helium abundance obtained 
by Koppen et al. I|1991|l is reportedly less accurate, due 
to incompleteness of the ionization correction. Regarding 
the heavy elements, the average discrepancy between our 



results and those by Koppen et al. I|1991f) is a factor 1.3, 
which means excellent agreement, for all elements and ob- 
jects well measured by them, that is, those not marked 
with a colon (:). The remaining elements, for which their 
abundances are reportedly less accurate, present a much 
larger deviation, so that we believe our measurements do 
represent an improvement in the knowledge of these ob- 
jects. 



In the second study by Cuisinier et al. (|1996(l . fea- 
turing nebulae located above the galactic plane, chemi- 
cal abundances have been derived for the objects Ml-9, 
Ml-16 and YC 2-5. Considering the first two, the He/H 
abundances agree with our values within 20%, and the 
heavy elements show an average discrepancy of a factor 
1.8, which is especially large for argon. For YC 2-5 the 
discrepancy is higher, about a factor 2.5, but it should 
be mentioned that Cuisinier et al. 1|1996H could not derive 
the sulfur abundance and their argon abundance is uncer- 
tain due to incompleteness of the ionization correction. On 
the other hand, we could not measure the nitrogen abun- 
dance, which is indeed very low according to Cuisinier et 



al. H1996J) . Also, its odd postion in Fig. 2 could be con- 
ciliated if its oxygen abundance is lower, as suggested by 
Cuisinier et al. I|19961) . Clearly, this object deserves further 
study. 

The comprehensive work of Kingsburgh & Barlow 
Q1994I hereafter KB94), based on optical spectrophoto- 
metric observations and UV data from the IUE satellite. 
Five nebulae from Table 1 (IC 2003, J900, Ml-8, Ml-13, 
and M3-3) are included in the KS94 sample, for which a 
generally good agreement is observed in the abundances 
within the uncertainties of each work. This is particularly 
interesting, as the abundances derived by KB94 use in- 
dependent ICFs, which are based on detailed photoion- 
ization models. For the He/H ratio, the abundances of 
these objects derived by KS94 show an average deviation 
of A(He/H) = 0.018 relative to our data, which corre- 
sponds roughly to 15%. The worst case is J900, for which 
a discrepancy of about 30% is obtained, but as we will see 
later, our value is in very good agreement with the recent 
result by Kwitter et al. J2003). For the heavy elements, 
the average deviations (dex) of the KS94 sample relative 
to our data are Ae(O) = 0.14, Ae(N) = 0.15, Ae(Ar) = 
0.27, and Ae(S) = 0.32, which reflects the fact that the 
sulfur and argon abundances are generally less well known 
than the oxygen and nitrogen data, in view of their lower 
relative abundances. The average discrepancy factors for 
oxygen and nitrogen are about 1.5, which are quite reason- 
able. For argon the discrepancy factor is somewhat larger, 
roughly 2.2. For sulfur we get a factor of 1.4 for the ob- 
jects IC 2003, Ml-8 and Ml-13. For J 900 and M3-3 the 
agreement of the results by KB94 and our data is poor, 
but the sulfur abundances of these objects were not well 
determined in the work of KB94. 

Perinotto & Corradi (1098.) analyzed the chemical 
structure of some bipolar planetary nebulae by means of 
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long slit spectrophotometry, considering several regions 
inside the nebulae in order to detect spatial abundance 
variations through the nebulae. There are three objects in 
common with our work, namely Ml-13, Ml-16 and M3-2. 
The He/H abundances for these objects are similar to our 
values within 20%, especially the relatively large abun- 
dance found in M3-2, which is close to our result and also 
to the data of Koppen et al. (J1991J1 . showing that this ob- 
ject is really helium-rich, as can be seen in Figure 1. The 
abundances of the heavy elements O, N, S and Ar are in 
agreement with our values within a factor of 2, in average, 
with a better agreement for Ml-16 and M3-2 and a poorer 
one for Ml-13, in which case our results are closer to those 
by KB94. 

Kholtygin (1998) determined oxygen abundances of 
about 70 planetary nebulae taking into account the ef- 
fect of small temperature and density fluctuations on the 
line intensities of the nebular lines. The objects IC 2003, 
IC 2165, J320 and J900 are included in his sample, for 
which the oxygen abundances are e(O) = 8.70, 8.50, 8.58 
and 8.45, respectively. The agreeement with our results 
is excellent, within 10% for the first three nebulae and of 
a factor 1.7 for J900. The case of IC 2165 is particularly 
interesting, as we have mentioned the discrepancy of the 
helium abundance and electron temperature relative to 
some previous work. 

The objects J320 and J900 were part of a recent survey 
of PN which aimed at providing a homogeneous spectro- 
scopic database in the wavelength range of 3600-9600 A 
(Milingo et al. l2UU21 Kwitter et al. l2UU3|) . This survey had 
the advantage of using the near infrared lines of [SIII] A 
9069, 9532 A, so that their abundances are expected to be 
more accurate, especially for sulfur, in view of a better de- 
termination of the ionic abundances of S +2 . A comparison 
of our results for J900 (Kwitter et al.|2003) shows a very 
good agreement for all elements better than 0.2 dex, or 
with an average discrepancy of a factor 1.3. This fact sup- 
ports our data for this nebula in view of the discrepancies 
found relative to some of the earlier work already quoted, 
especially for nitrogen (Aller et al. I1983|) . argon (Koppen 
et al. I1991f) . oxygen fKholtvgin I1998|) . and the elements 
argon and sulfur (Kingsburgh & Barlow l*L994|l . For J320 
(Milingo et al. l2002|> the agreement is good for oxygen and 
argon, with an average discrepancy of a factor 1.3, that is, 
the abundance differences are within 0.2 dex. For sulfur, 
we find an abundance higher by a factor 2, or Ae ~ 0.30 
dex. In fact, both our results for the sulfur abundances 
of J320 and J900 confirm the remarks by Kwitter et al. 
(2003) and Milingo et al. (|2002fl . in the sense that use of 
the optical lines only probably leads to an overestimate of 
the sulfur abundance. For nitrogen, there is no agreement 
as this element is probably underestimated in our sample. 

4. The radial abundance gradient 

One of the most interesting consequences of the determi- 
nation of abundances of PN located in the galactic an- 
ticenter is the effect they might have on the determina- 



tion of the space variations in the radial abundance gra- 
dients. An analysis by Maciel & Quireza (|19991 hereafter 
MQ99) suggested that the gradients probably flatten out 
at galactocentric distances larger than R ~ 10 kpc (assum- 
ing i? = 7.6 for the galactocentric distance of the LSR), 
especially for the O/H and Ne/H ratios. It is therefore 
interesting to investigate whether the presently derived 
abundances follow the same trend. 

We have then produced a merged sample including the 
objects of MQ99 and the present results. Both samples 
form a homogeneous group, as the abundances in MQ99 
were obtained using basically the same methods as here. 
In fact, all abundances are generally similar within the 
average uncertainty of about 0.1 to 0.2 dex. 

In Figure 01 we plot our new results for oxygen (open 
circles) along with the nebulae presented by MQ99 (filled 
circles). We have not included the object K3-70, for which 
apparently no reliable distance has been determined. We 
adopted the same galactocentric distance of the LSR as 
MQ99 and the PN distances are from Cahn et al. (|1992|l . 
except for H3-75, for which we have used the distance 
given by Amnuel et al. 1|1984|> . The adopted heliocentric 
distances d and the galactocentric distances R (kpc) are 
shown in the last two columns of Table 1 . Use of a differ- 
ent distance scale such as the catalogue by Maciel (1984 ) 
produces essentially the same results. As has been shown 
before (see for example Maciel & Koppen 119941 hereafter 
MK94), the use of different statistical distances does not 
significantly affect the magnitudes of the gradients. 

The results indicate that the radial abundance gradi- 
ent, which is reasonably constant between galactocentric 
distances from 4 to 10 kpc, decreases for greater distances. 
The effect of the anticenter nebulae of the present work in 
the determination of the gradient is clearly seen in Figure 
0] as our sample was chosen specifically to cover galacto- 
centric distances greater than Rq, so they fill the lack of 
data for these distances in previous samples, such as MQ99 
and MK94. This can also be seen by the least squares fits 
shown, both for a linear correlation (the straight line in 
Figure 0J and a quadratic regression (curve). For the lat- 
ter, the inner gradient is about —0.09 dex/kpc for galacto- 
centric radii between 4 and 5 kpc, reaching an essentially 
flat gradient for R > 11 kpc. The average slope as in- 
dicated by the straight line is —0.05 dex/kpc, which is 
slightly lower than previous determinations for galacto- 
centric distances close to Rq. The flattening at large dis- 
tances is seen by the difference between the straight line 
and the parabola, especially at R < 5 kpc and R > 12 kpc. 
This result supports other observational results for HII re- 
gions (see for example Vflchez & Esteban fT996) an< ^ some 
theoretical models for the evolution of spiral galaxies (see 
Pagel l*P9"9"7l and references therein). 

It should be noted that the dispersion of the abun- 
dances increases at large galactocentric distances, proba- 
bly due to the fact that the nebulae are located farther 
away than solar neighbourhood objetcs, which are better 
studied. For the remaining elements a similar behaviour is 
observed, which is also seen in comparison with the earlier 
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Fig. 4. The O/H gradient from the merged sample of Maciel & Quireza {1999, filled circles) and the present sample 
(open circles). Linear and quadratic least squares fits are also shown. 



sample of Maciel & Koppen ( 1 !)!) 1 }. For a detailed discus- 
sion of these effects the reader is referred to the recent 
work by Maciel & Costa l|2U03|l and Maciel et al. PUU5|) . 

It is interesting to compare the variation of the O/H 
gradient with galactocentric radius as indicated by Figure 
01 with some recent results derived from different galac- 
tic objects. Data from O, B stars as discussed by Smartt 
H2000I see also Rolleston et al. 120001 suggest a constant 
slope for all galactocentric radii. The different behaviour 
of both samples is probably related to the time variation 
of the O/H gradient, as recently discussed by Maciel et 
al. I|2003[) . According to that work, the average gradient, 
which can be roughly defined as the gradient from galac- 
tocentric radii of 4 kpc to about 10 kpc, has flattened out 
from —0.08 dex/kpc to about —0.06 dex/kpc in the last 
5 Gyr, allowing for an average uncertainty of about 0.2 
dex/kpc. Therefore, it is tempting to conclude that at the 
time when the oldest disk PN have been formed - sev- 
eral Gyr in the past - the average gradient was somewhat 
steeper than indicated by younger objects such as O, B 
stars and HII regions. However, some flattening was al- 
ready present near the outer edge of the Galaxy, probably 
reflecting the lower star formation rate in that region as 
compared with the regions closer to the galactic center. 
Such flattening is apparent from our PN data, as can be 
seen in Figure As the galactic evolution proceeded, the 
flattening initially located near the outer Galaxy spread 
towards the more central regions, which explains the lower 
and approximately constant gradient shown by the O, B 
stars. Since the temporal flattening rate has decreased in 
the last few Gyr, as suggested by the results of Maciel et 
al. H2003f) . the gradients presented by young objects such 
as O, B stars and HII regions are expected to be similar, 
which is supported by the observations, again within an 



average uncertainty of 0.02 dex/kpc. Such a scenario is 
also supported by recent work on HII regions, for which 
gradients as low as about —0.04 dex/kpc have been pro- 
posed (Deharveng et al. [2000). On the theoretical side, 
gradients that flatten out in time with a decreasing flat- 
tening rate in the last few Gyr are supported by the recent 
models proposed by Hou et al. (2000). 
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